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Seismic  Hazard  Study  for  Utah 


I.  INTRODUCTION 


Several  of  the  proposed  deployment  modes  for  the  MX  missile  system  have  in¬ 
volved  basing  of  the  system  in  the  alluvial  basins  of  Nevada  and  western  Utah.  In 
support  of  the  MX  missile  program,  the  Terrestrial  Sciences  Division  of  the  Air 
Force  Geophysics  Laboratory  has  undertaken  a  program  to  evaluate  the  seismic 
hazard  with  these  two  states.  In  this  report,  results  of  the  analysis  for  the  Utah 
area  are  presented.  Doth  probabilistic  and  deterministic  approaches  to  the  esti¬ 
mation  of  seismic  hazard  have  been  employed.  The  results  of  the  study  are  set 
forth  as  contour  plots  of  maximum  credible  and  probabilistic  peak  ground  motion 
estimates  for  hard  rock  sites.  In  addition,  annual  risk  curves  are  presented  for 
risk  assessments  in  Utah.  No  attempt  has  been  made  to  modify  the  rock  site 
accelerations  to  compensate  lor  the  deep  alluvial  cover  typical  of  the  potential  MX 
basing  sites. 

The  approach  taken  in  this  study  was  to  produce  a  conservative,  though  realis¬ 
tic,  estimate  of  the  seismic  hazard  in  the  area  of  interest.  However,  the  long-term 
seismicity  patterns  in  this  region  are  not  well-defined  and  subject  to  debate.  To 
provide  some  sense  of  the  uncertainty  associated  with  the  analysis,  two  distinct 
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models  of  the  regional  seismicity  have  been  used  in  the  probaoilistic  evaluation. 
The  results  from  both  assumptions- are  presented  and  compared. 


2.  UTAH  SEISMICITY 

The  earliest  record  of  seismic  activity  in  Utah  is  of  a  tremor  felt  in  central 

Utah  in  1853.  *  Since  that  event,  earthquakes  of  sufficient  size  to  be  felt  have 

2 

occurred  within  the  state  at  the  rate  of  approximately  one  per  year.  During  the 
historical  era,  earthquake  epicenters  in  Utah  have  been  concentrated  in  a  distinct 
band,  running  from  the  southwestern  corner  of  the  state  to  north-central  Utah 
(Figure  1).  Although  the  preinstrumental  earthquake  catalog  is  biased  in  favor  of 
this  zone  by  virtue  of  the  population  distribution  during  most  of  the  state's  history, 
instrumental  coverage  in  the  last  two  decades  has  re-emphasized  this  pattern  of 
seismic  activity.  However,  significant  earthquakes  have  also  occurred  widely 
throughout  the  state  and  outside  of  this  trend. 

The  largest  earthquake  within  the  state  was  the  1934  Hansel  Valley  earthquake 
that  had  a  reported  magnitude  of  M=6.  6.  This  event  was  centered  just  north  of  the 
Great  Salt  Lake.  At  least  57  earthquakes  have  been  recorded  in  Utah  between  1853 
and  1978,  causing  at  least  minor  damage.  Modified  Mercalli  intensity  V  or  greater. 
With  use  of  the  earthquake  catalog  for  all  of  Utah  from  1850  to  1974,  evaluation  of 
recurrence  curves  for  Utah  gives  an  incremental  curve  of 

log  N  =  1.  15  -  0.  76  M  (1) 


and  a  cumulative  curve  of 


log  N  =  1.47  -  0.  79  M  (2) 

2 

where  N  is  the  number  of  events  per  year  per  1000  km  of  magnitude,  M,  for  the 

,  *  4 

met  ..mental  curve  or  equal  to  or  greater  than  M  for  the  cumulative  curve. 

1.  Cook,  K,  L, ,  and  Smith,  R.  B,  (1967)  Seismicity  in  Utah,  1850  through  ,lune 

1975,  Bull.  Seismol.  Soc.  Am  57:689-718. 

2.  Cook,  K.  L.  (1972)  Earthquakes  along  the  Wasatch  Front,  Utah  -  The  record 

and  the  outlook,  in  Environmental  Geology  of  the  Wasatch  Front,  1971, 

L.  S.  Hilpert,  ed. ,  Utah  Cool.  Assoc.  ,  Publ,  1,  pHl-1129. 

3.  Meyers,  H. ,  and  von  Hake,  C.  A.  (1976)  Earthquake  Data  File  Summary. 

National  Geophysical  and  Solar -Terrestrial  DataCenter,  Report  KCf\D-5, 

4.  United  States  Geological  Survey  (1976)  A  Study  of  Earthquake  Losses  in  the 

Salt  Lake  City,  Utah  Area,  U.  S.  Geological  Surv,  Open  File  Report  76-89. 
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Earthquakes  in  the  Utah  area  tend  to  be  shallow  events,  with  the  majority  of 

5 

epicenters  having  depths  above  10  km.  However,  only  the  1934  Hansel  \  alley 
earthquake  can  be  associated  with  surface  faulting.  An  earthquake  swarm  in 
southwestern  l  tali  in  1971  produced  fractures  in  alluvium  that  might  be  related  to 
faulting.  '’  In  general,  the  pattern  and  locations  of  earthquake  epicenters  in  l  tab 
tend  to  be  diffuse  and  difficult  to  associate  with  specific  faults.  ^  Although  faults 
are  common  throughout  most  of  Utah  (Figure  2),  the  age  of  most  recent  faulting 
is  unknown  on  many  of  them.  ‘  In  the  assessment  of  earthquake  hazard,  faults 
having  Quaternary  displacements,  those  younger  than  2  million  years,  are  often 
considered  potentially  active.  Using  the  limited  data  available,  Anderson  and 
Miller  have  compiled  a  map  of  Quaternary  faults  in  Utah  (Figure  3).  They  consider 
this  map  to  be  preliminary;  additional  faults  are  likely  to  be-  reported  with  further 
geologic  studies.  The  distribution  of  known  Quaternary  faulting  supports  the 
general  conclusion  attained  from  the  historical  earthquake  record:  That  while  most 
seismic  activity  is  associated  with  the  Wasatch,  Eisinore,  and  Hurricane  Fault 
zones,  potc.  tial  exists  for  significant  seismic  activity  in  other  areas  of  the  state. 
Historically,  the  1934  Hansel  Valley  earthquake  marks  the  largest  event  oc¬ 
curring  in  Utah.  However,  it  is  possible  that  larger  earthquakes  will  occur  there 
in  the  future.  Detailed  studies  of  the  geologic  evidence  of  Holocene  fault  displace¬ 
ments  along  the  Wasatch  Front  suggest  that  individual  surface  faulting  events  have 

9 

produced  displacements  in  the  range  of  0.  8  to  3.  7  m.  Given  the  fault  displacement, 
D,  the  causative  event  magnitude,  M,  can  be  estimated  by  the  functional  relationship 
derived  for  normal  faults.  10  The  Wasatch  Fault  zone  can  then  be  assumed  to  have 


5.  Smith,  R.  H. ,  and  Sbar,  M.  L.  (1974)  Contemporary  tectonics  and  seismicity 
of  the  Western  United  States  with  emphasis  on  the  inter -mountain  seismic 
belt,  Gcol,  Soc,  ■■Urn.  Hull.  72:1205-1218. 

G.  Hosier,  D,  I,,  and  Smith,  R.  B.  (1982)  Seismic  moment  rates  in  the  Utah 
region.  Bull.  Seismol.  Soc.  Am.  72:525-551. 

7,  Anderson,  L.  W. ,  and  Miller,  D.  G.  (1979)  Quaternary  faulting  in  Utah,  in 

Earthquake  Hazard  along  the  Wasatch  and  Sierra  Nevada  Frontal  Fault 
Zones,  U.  S.  Geol.  Surv.  ,  Open  File  Heport  80-801,  pp.  1 94 -226, 

8.  Allen,  C.  R,  (1975)  Geologic  criteria  for  evaluating  seismicity,  Geol.  Soc. 

Am,  Bull.  86:1041-1057.  ~ 


9.  Swan,  F.  H, ,  Schwartz,  D.  F. ,  and  Cluff,  L,  S,  (1979)  Recurrence  of  surface 
faulting  and  moderate  to  large  earthquakes  on  the  Wasatch  Fault  zone  at 
Kaysville  and  Hobble  Creek  sites,  Utah,  in  Earthquake  Hazards  along  the 
Wasatch  and  Sierra  Nevada  Frontal  F  ault  Zones.  U.  S.  Geol.  Surv.,  Open 
File  Report  60-801,  pp.  227-275. 

10.  Slcmmons,  D.  B.  (1977)  Faults  and  Earthquake  Magnitude,  U.  S.  Army  Engin. 
Waterways  Experiment  Station,  Miscellaneeus  Paper  S-73-1. 
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Figure  2.  Known  and  -Suspected  Faults  in  l  tab.  Known  faults,  solid  lines; 
suspected  faults,  dashed  lines.  "After  Cook  and  Smith  (1967). 


supported  earthqu  ikes  of  magnitude  M  =  (>.  7  to  7.1  repeatedly  during  the  Holocene, 
approximately  the  last  10,000  years.  It  must  he  assumed  that  an  earthquake  of 
about  magnitude  \  17 .  a  could  occur  in  l  tali,  at  least  along  the  Wasatch  Fault  .rone. 
This  value  is  commonlv  taken  as  the  maximum  credible  earthquake  for  Ft  ah. 

The  historical  record  of  earthquakes  implies,  through  Fq.  (2),  the  recurrence 
of  a  major  earthquake  of  magnitude  M  =  7.0  or  greater,  approximately  once  every 
52  years  in  Utah,  or  130  years  for  an  M=7.  5  event.  Based  on  several  assumptions, 
as  well  as  extrapolation,  the  geologic  evidence  suggests  that  earthquakes  of  this 


1  1 


order  of  magnitude  should  have  return  periods  between  50  and  400  vears  along  the 

t  .  ‘i 

Wasatch  Fault  zone.  During  the  last  LiG  vears,  onlv  the  Hansel  \  allcv  event 
approached  this  level  of  magnitude.  In  geologic  terms,  the  historical  record  is 
extremelv  short  and  extrapolation  Irom  these  data  tor  magnitudes  which  have  not 
been  recorded  is  probabilistic.  Analysis  of  earthquake  catalogs  for  China,  .lapan, 
and  the  Middle  East,  covering  periods  between  2,000  and  1,000  years,  indicate 
that  spatial  and  temporal  variations  in  seismicity  of  periods  equal  to  or  longer  than 
the  Utah  historical  record  are  common.  8|  1  sing  various  windows  of  the  I  tali 

catalog,  Arabasz  et  al  have  shown  that  return  periods  ol  between  .4  and  1,  405  vears 
could  be  estimated  for  a  7.  5M  earthquake  on  the  Wasatch  front,  Tvpic  tllv,  in 
estimate  of  the  order  of  several  hundred  vears  is  assumed  for  the  maximum  cred¬ 
ible  earthquake  of  M  =  7.  5. 

:s.  SK1SM1C  KF.EION  AI.IZATION 

Subdivision  of  the  t  tali  area  into  seismic  source  regions  is  required  in  order 
to  analyze  the  seismic  hazard  Ln  the  state.  The  goal  of  this  regionalization  is  to 
produce  a  set  of  seismic  source  regions  whose  seismicity  is  homogeneous  and 
well-defined.  Limitations  imposed  bv  short  seismic  history,  the  diffuse  nature  of 
the  seismicity  in  the  area  of  study,  and  questions  concerning  the  tectonic  procc  =ses 
of  the  region  --  all  affect  the  process  of  defining  these  source  regions.  In  this 
section,  two  regionalizations  are  made  in  an  attempt  to  compensate  tor  the  uncer¬ 
tainties  previously  stated. 

The  n*  tin  trend  of  seismic  activity  through  the  state  is  associated  •  it li  ’lie 
Intermountain  Seismic  Holt  (ISB),  a  trend  of  earthquakes  ranging  from  northwestern 
Arizona  north  to  the  Montana-British  Columbia  border  (Figure  41.  ’  This  zon<  is 
one  of  the  most  active  seismic  regions  in  the  continental  l  iiitcd  States.  It  has  been 
suggested  that  the  [SB  defines  the  active  margins  of  several  sttbplate.u  of  dt"  North 
American  Plate.  Hie  suhplate  margins  are  believed  to  correlate  csscntiallv  with 
the  physiographic  province  boundaries  that  arc  also  shown  in  Figure  4.  Srisetiritv 

11.  McGuire,  H,  K.  (1S'7!'1  Adequacy  of  simple  proltabilitv  models  lor  calculating 

felt  shaking  hazard,  using  the  Chinese  earthquake  catalog,  Hull,  Sejs mol, 
Soc,  Am,  f)‘ld!77-:i!12. 

12.  York,  .1.  E.  ,  Cardwell,  I?.,  anti  M,  J.  <  lt'7()l  Seismieit  v  and  Quatcrnam 

faulting  in  China,  Bull,  Seismol,  Soc,  Am.  (>6:lt'!H-2002. 

F!.  Arabasz,  W.  .1.,  Smith,  It.  B. ,  an*  I  Hie  bins,  W,  I).  (lhT’h  Earthquake  studies 
along  tin'  Wasatch  Front,  l  tali:  network  monitoring,  seismicity,  and  sei-mic 
hazards,  in  Earthquake  Hazards  along  the  Wasatch  ami  Sierra  Nevada 
Frontal  Fault  Zones,  F.  S.  Geol.  Sun.,  Open  File  He  port  B0  —HOI .  pp.  l-:!.l. 
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Ki L'(J r<*  •}.  inr<*  rmountain  S I'iyivic  Molt  (stipplod  ar*n)  and  I ’  \  >10,' :  iphii' 
1‘rovinros  of  the  I  t.d>  \i«  .  mm  \ft«r  and  Sh.u  (U'7-l), 


•  s  itt-.in  tin'  ISM  >pp.  u's  to  result  from  lorn'-  associated  with  the  interaction  of  the 
North.  American  uvl  Pacific  Plates  ami  pos.-inlv  modified  by  a  hypothesized  mantle 
plume  prescntlv  centered  at  Yellowstone.  In  the  L'tah  area,  the  ISM  divides  the 
Basin  .inti  Range  province  on  tin  west  from  the  Middle  Kocky  Mountains  and  the 
Color  ado  Ml ateau  on  tin-  east.  The  distinct  tectonics  of  each  of  these  provinces  is 
siunifietuu  in  tin-  seismie  regionalization  of  the  area. 

Kpcirogenic  uplift  has  occurred  in  all  three  provinces,  but  it  is  the  dominant 
tectonic  force  onh  in  the  Colorado  Plateau  and  the  Middle  llockv  Mountains.  In 
tlte  Colorado  Plate  tu,  the  uplift  has  been  uniform  throughout  the  province,  with 
l i 1 1 1 •  crustal  deformation.  Tile  rough  topography  of  the  plateau  is  primarily  the 
result  of  fluvial  >  rosion  of  undisturbed  sedimentarv  deposits.  The  widespread 
ni'.ure  ol  tin-  uplift  and  lack  of  crustal  deformation  suggest  a  relatively  stable  zone 
ot  low -R  .cl  s.  ismicif.  th  tt  would  be  uniformly  distributed;  in  fact,  concentrations 
of  seismic  microactivity  that  have  occurred  in  this  region  appear  to  be  associated 
with  human  activity  such  as  minim..  In  the  Middle  Rocky  .Mountains,  the  epeiro- 
gvnie  uplift  has  been  accomplished  hv  extensive  normal  faulting  and  mountain  forma¬ 
tion.  (in  the  scale  of  interest  in  this  report,  the  seismic  activity  would  also  be 
expected  to  be  uniformly  distributed,  but  the  region  would  be  expected  to  have  a 
higher  level  of  seismicity  than  the  Colorado  Platorut.  Each  of  these  provinces  can 
be  considered  a  single  source  region  with  uniform  potential  for  seismic  activity. 

Vn  exception  is  the  margin  of  the  Colorado  Plateau  abutting  the  liasin  and  Range; 
this  will  be  discussed  later. 

1  he  major  tectonic  process  in  the  Basin  and  Range  is  Unilateral  crustal  exten¬ 
sion  which  has  been  occurring  for  the  last  30  million  years.  The  rifting  has 
continued  episodically  to  the  present  day.  although  the  orientation  and  center  of 
active  extension  has  changed.  During  this  period,  tvpiral  estimates  of  the  amount 
of  spreading  range  from  f>0  to  .iOO  km,  with  the  most  common  estimate  being  about 
100  km.  '  The  spreading  has  been  accomplished  by  extensive  normal  faulting, 
graberi  formation,  and  voleanism.  Uthough  probably  related  to  the  broad  uplift 
of  the  region,  the  exact  relationship  between  uplift  and  extension  of  the  Basin  and 
Range,  as  cause  or  effect,  is  uncertain.  *  ’ 

It.  Smith,  R.  B. ,  Winkler,  P.  I...  Anderson,  J.  G.  ,  and  So  liol  z  (1974)  Source 
mechanisms  of  micro-earthquakes  associated  with  underground  mines  in 
eastern  l'tah.  Bull.  Scismol.  Soc.  Am,  <>4:1295-1317. 

15.  Baton,  G.  P.  (1970)  A  plate  tectonic  model  for  late  Cenozoic  spreading  in 

the  Western  l  nited  States, 

R.  K.  Riecker,  ed.  ,  Am, 

16.  Thompson,  G.  A.,  and  Bruce,  1).  B.  (1974)  Regional  geophysics  of  the  Basin 

anil  Range  province,  in  Annual  Review  of  Barth  and  Planetary  Sciences. 

T.  A.  Donath,  oil..  Annual  Reviews,  Inc,,  Palo  Alto,  California,  pp.  213-238. 


in  itio  uranue  itxrt;  tectonics  anti  lViagmans: 
Gcophys.  I  nion,  Washington,  D.  C. ,  pp.  7-31 
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A  question  with  significant  implications  for  seismic  regionalization  is  the 
distribution  of  active  extension  in  the  Basin  and  Range.  During  historic  times, 
most  significant  earthquakes  in  the  province  have  occurred  in  very  limited  zones 

that  tend  to  he  at  the  margins  of  the  province.  Late  Quaternarv  faulting  in  the 

17 

region  is,  however,  very  uniformly  distributed  as  is  low  level  seismicity.  At 
least  three  hypotheses  have  been  advanced  in  an  effort  to  explain  this  apparent 
discrepancy,  first,  the  historical  record  is  a  short-term  distribution  of  earth¬ 
quakes  in  the  province;  second,  the  seismic  activity  (and  presumably,  extensional 
activity)  is  concentrated  in  very  limited  zones  at  any  one  time;  these  zones  migrate 
throughout  the  Basin  and  Range.  *  Return  periods  for  large  earthquakes  in  any 
specific  area  are  thought  to  be  of  the  order  of  thousands  of  years.  Finally,  the 
last  c  oncept  assumes  that  extensional,  as  well  as  seismic  activity,  is  confined  to 

the  Basin  and  Range  margins  :ind,  on  the  eastern  edge,  could  be  gradually  migra- 

5 

ting  into  the  Middle  Rocky  Mountains  and  Colorado  Plateau. 

Historical,  as  well  as  geologic  evidence,  is  insufficient  to  determine  the 
correct  analysis  of  tin'  seismicity  of  the  Basin  and  Range,  Each  hypothesis,  how¬ 
ever,  would  require  a  different  approach  to  the  seismic  regionalization  of  the 
province.  The  first  concept  would  imply  that  the  Basin  and  Range  could  be  con¬ 
sidered  as  a  single  or  at  least  a  limited  number  of  zones  having  uniform  seismic 
activity.  I  sing  the  second  hypothesis,  one  would  have  to  delineate  any  area  with 
indications  of  recent  major  earthquakes  as  a  potential  earthquake  source.  As  the 
future  migration  of  the  active  zones  cannot  be  predicted,  the  remaining  area  would 
be  assumed  to  have  a  uniform  probability  of  earthquake  occurrence.  On  the  scale 
of  regionalization  used  in  this  study,  this  would  degenerate  into  the  same  conditions 
as  the  first  approach.  On  the  basis  of  the  last  thesis,  seismic  activity  would  be 
confined  to  the  province  margins  with  low -level  background  seismicity  in  the 
interior  of  the  Basin  and  Range'. 

It  is  concluded  that  two  models  of  seismic  regionalization  will  provide  an 
adequate  range  of  seismic  hazard  estimations  for  the  t  tah  area.  In  the  first, 
known  as  the  t  niform  Seismicity  (LS)  model,  a  major  earthquake  is  assumed  to 
have  uniform  probability  of  occurrence  within  several  gross  regions;  in  the  second, 
the  Subplate  Margin  Seismicity  (SMS)  model,  the  major  seismic  source  regions  are 
assumed  to  align  with  the  historical  seismicity  patterns  of  the  region.  The  primary 
distinction  is  in  the  treatment  of  the  Basin  and  Range  province  earthquake  potential. 


17,  Slemmons,  I).  B.  (H’67)  Pliocene  and  Quaternary  crustal  movements  in  the 

Basin  and  Range  Province,  I  SA,  .1  our nal  Geosciences,  Osaka  City  Univ.  , 

10:  m  -  10B.  ‘ 

18,  Ryall,  A.  (1077)  Earthquake  hazard  in  the  Nevada  region.  Bull.  Seismol.  Soe. 

Am,  67:517-5. id. 
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Figure  5.  I  niform  Seismic itv  Model  Regionalization  [(Modified  From  Gieensfelder 
et  al.  (1080)] 


seismicity  parameters  are  based  primarily  on  the  work  of  Greensfelder  ct  al. 
Additional  analysis  of  seismicity  was  conducted  for  the  Colorado  Plateau-Rocky 
Mountains  region  to  the  east  of  the  area  included  in  that  study.  The  regionalization 
is  based  on  late  Cenozoic  patterns  of  deformation,  relative  Holocene  strain  rates, 
and  instrumental  historic  seismicity.  Greensfelder  et  al,  conclude  that  Holocene 
strain  rates  and  historic  seismicity  are  in  rough  agreement  with  differences  of  less 
than  a  factor  of  :i.  Other  investigators  have  suggested  that  the  differences  are 
primarily  due  to  the  methods  used  for  estimating  strain  rate  rather  than  changes 
in  long-term  seismicity.^ 

Based  on  post -instrumental  seismicity,  1!T'12  to  lf'T.'l,  the  seismicity  param¬ 
eters  for  each  region  were  calculated  (See  Table  1).  The  values  given  for  regions 
1  through  6  are  based  on  earthquakes  west  of  region  6  to  10HJW.  1’he  method  of 
analysis  for  region  7  was  similar  to  that  used  for  the  other  regions. 


Table  1.  Triform  Seismicity  Model  Parameters 


*N  -  number  of  events  of  maqnitudc  M  or  qreater  per  year  let  km'  . 


Greensfelder,  11.  W. ,  Kintzer,  T.  ,  and  Somervill,  M.  K.  (1MJ10)  Seismo- 
tectonic  regionalization  of  the  Great  Basin,  :tnd  comparison  of  moment  rat 
computed  from  Holocene  strain  and  historic  seismicity:  summary,  Gcol. 
Soc,  Am.  Bull,  <H:r>lil-f>2S, 


tn  addition  to  recurrence  curve  information,  maximum  magnitude  earthquakes 
are  also  assigned  for  each  region  and  shown  in  1'able  1,  These  values  are  based 
on  interpretation  of  geologic  evidence  or  the  largest  event  to  occur  in  historical 
times.  In  general,  modifying  this  value  has  a  significant  effect  only  on  the  hazard 
estimation  for  very  long  return  periods.  Unless  the  maximum  magnitude  earth¬ 
quake  is  grossly  in  error,  probabilistic  risk  estimates  in  the  range  of  interest  for 
this  report  will  not  be  modified  substantially  by  changing  this  parameter. 

3.2  Suhplate  Margin  Seismicity  Model 

In  this  model,  the  regional  seismicity  is  assumed  to  concentrate  at  the  suh- 

r, 

plate  margins,  as  suggested  by  historical  seismic  activity.  The  basic  regional¬ 
ization  of  the  uniform  seismicity  model  is  maintained  outside  of  Utah,  The  major 
distinction  between  this  model  and  the  previous  one  is  the  lowering  of  the  seismic 
potential  in  region  3  :uui  the  redefining  of  region  6  bv  subdivision  into  two  new 
zones,  with  alteration  of  the  activity  levels  in  each  new  zone  (figure  6).  The  new 


figure  6.  Subplate  Margin  Seismieitv  Moiled  Regionalization 

10 


region  t>l)  corresponds  to  the  Wasatch  Front  Study  Area  defined  by  Arab  as  z  et  al. ; 
it  is  an  areal  expansion  of  the  equivalent  segment  of  region  6  from  the  uniform 
■seismicity  model.  '  Adjustments  in  the  abutting  regions  3  and  7  are  required 
by  this  change. 

The  recurrence  function  used  for  the  tu  w  region  a  was  derived  from  a  previous 

2  0 

study  lor  :ut  area  of  Nevada  closely  resembling  the  redefined  region  3.  It  is  felt 
that  the  'ifferera  area  used  in  this  evaluation  w  ould  not  introduce  errors  of  a  magni¬ 
tude  approaching  those  expected  from  the  statistically  short  duration  of  the  earth¬ 
quake  sample.  For  regions  6a  aid  6i>,  speeifie  analysis  of  seismic  recurrence  for 
each  was  conducted  using  the  instrurm  atal  earthquake  catalog  from  1932  to  1078. 
This  included  an  analysis  of  eomoleter.ess  of  the  catalog  to  ensure  st  Uistical 

*>  i 

stability.  In  addition,  the  modeled  maximum  m  ignitude  e:u  thquaKcs  in  regions  > 
uid  6a  were  reduced  in  order  to  proviiie  i  lower  level  of  seismic  activity  in  those 
region-  us  oompiu’ed  to  region  61  as  indicated  b\  the  historic  :1  record.  The  recur - 
re-nee  parameters  for  the  subpltte  margin  eismicit\  model  are  given  in  Table  2, 


fal.'le  2.  Suhplate  Margin  Seij.it  ieit;  Model  Parameters 
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20.  P.attis,  J.  C. ,  and  Hill,  K.  T.  (1977)  /Analysis  of  Seismicity  and  Tectonics 

of  the  Central  and  Western  United  States,  Texas  Instruments,  Inc.  .  Inter i m 
Scientific  Report  ALEX  (02) -iSR -7 7-01,  Dallas,  Texas. 


21.  Stepp,  J.  C.  (1972)  Analysis  of  completeness  of  the  earthquake  sample  in  the 

Puget  Sound  area  and  its  effect  on  statistical  estimates  of  earthquake  hazard, 
in  Proceedings  of  the  International  Conference  on  Microzonation  for  Safer 
Construction  Research  and  Application,  I'nlv.  of  Washington,  Seattle, 


1'ho  probabilistic  hazard  estimation  process  is  based  on  a  method  proposed  by 
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Cornell  and  implemented  in  a  1-' OUT  HAN  computer  program  by  McGuire. 
l  sing  this  program,  we  found  that  the  temporal  and  spatial  distribution  of  seismic 
activity  is  combined  into  a  single  statement  of  the  probability  of  reaching  or  exceed¬ 
ing  a  given  ground  motion  level  at  any  site  of  interest.  In  addition  to  the  statistical 
models  of  seismicity,  given  in  the  preceding  section,  knowledge  of  ground  motion 
attenuation  in  the  region  of  interest  is  required.  In  the  following  paragraphs,  the 
statistical  method,  ground  motion  attenuation,  and  results  of  the  analysis  for  Utah 
arc  discussed. 

t.l  Risk  Kshm Jt mil  Method 

In  the  procedure  used  for  this  study,  earthquake  occurrence  within  each  sub¬ 
division  is  considered  a  Poisson  process.  This  implies  that  the  occurrence,  in 
space  and  time  within  any  subregion,  is  independent  of  any  preceding  event.  There 

is  strong  evidence  that  this  is  not  true  and  that  future  earthquake  occurrence  is  con- 

1 8  24 

neeted  to  past  seismic  activity  in  the  region.  *  “  However,  the  use  of  the 
Uoissonian  assumption  is  conservative;  it  does  not  introduce  additional  restrictions 
which  at  best  could  only  be  based  on  extremely  limited  data  samples. 

I  sing  the  Poisson  process,  we  note  for  each  source  region  the  probability 
that  the  ground  motion  will  reach  or  exceed  a  specified  level,  a^,  defined  as  the 
integral  of  the  product  of  three  parameters:  the  independent  probability  density 
functions  for  occurrence  of  an  earthquake  of  magnitude  S,  f<.;  the  density  function 
that  it  will  occur  at  distance  ft,  f^ ;  and  the  conditional  probability  that  given  an 
event  of  magnitude  s  and  at  distance  r,  the  ground  motion  will  reach  or  exceed  a^. 
This  can  be  stated  as: 

Pi  A  >  a  j-  ,  s'  Ft  A  >  a  I  s  and  rjf_  (s)f„(r)dsdrj  (2) 

g-  g  g-  g  s  K 

I  2:> 

where  ft.  A  >  a  |s  and  rj  is  the  conditional  probability." 


22.  Cornell,  C.  A.  (1968)  Engineering  seismic  risk  analysis,  Hull.  Seismol. 

Soc.  Am.  50:1503-1606.  . 

23.  McGuire,  K.  K.  (1976)  KOHTKAN  Computer  Program  for  Seismic  Risk 

Analysis,  P.  S.  Geol.  Survey,  Open  rile  Report,  ?6-67. 

24.  Sykes,  L.  K.  (1971)  /Aftershock  zones  of  great  earthquakes,  seismicity  gaps, 

and  earthquake  prediction  for  Alaska  and  the  Aleutians,  J.  Geophvs.  Kes. 
76:8081-8041. 
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The  function  fg(s)  is  derived  from  the  source  region  recurrence  curve  while 

f^(r)  incorporates  the  spatial  relationship  between  the  total  source  region  and  the 

site  of  interest.  The  conditional  probability  is  derived  from  a  ground  motion 

attenuation  function  discussed  in  the  next  section.  Evaluation  of  the  integral  yields 

the  probability  that  one  event  from  the  gi\  en  source  region  will  reach  or  exceed  a^. 

By  multiplying  this  value  by  the  expected  number  of  events  in  the  region  for  one 

year  and  summing  over  all  source  regions  with  some  predetermined  distance,  we 

obtain  the  total  expected  number  of  events  per  year,  or  annual  risk,  R^. 

The  risk  can  be  stated  in  alternate  forms;  for  example,  return  period,  P  ,  or 

lifetime  risks,  R...  The  return  period  is  given  bv  the  inverse  of  R  ,  and  corre- 
IN  A 

sponds  to  the  average  interval  of  time  between  the  site  of  interest  experiencing  the 
specified  acceleration  or  greater.  The  lifetime  risk,  R^,  is  the  probability  of  ex¬ 
ceeding  the  specified  acceleration  in  any  N  year  period.  1  or  a  Poisson  process, 
this  value  is  given  by: 

rn  =  1  -  ( 1  -R  a>N-  (4> 

As  an  example,  a  475-vear  return  period  ground  motion,  R^~  0.  0021,  is 
equivalent  to  R^q  =  0.  1.  This  is  equivalent  to  stating  that  in  any  50 -year  period 
the  level  of  ground  motion  having  a  475 -year  return  period  will  not  be  exceeded  at 
the  90  percent  confidence  level. 

4.2  Ground  Motion  Attenuation 

The  conditional  probability  density  function  discussed  in  the  preceding  section 
is  derived  from  ground  motion  attenuation  functions.  Various  empirical  studies 
have  been  conducted  to  derive  these  functions  (see  Ref.  23).  Typically  these 
equations  have  the  form 

g  =  a j  e‘V(R  +  as>'a4  (4) 

where  g  is  the  predicted  motion  level,  M  is  earthquake  magnitude  and  H  the  event- 
to-site  distance.  At  best,  these  equations  are  of  limited  accuracy  due  to  lack  of 
radiation  pattern  or  local  travel  path  and  site  geology  corrections,  and  the  standard 
deviations  associated  with  them  are  typically  large.  Most  equations  are  hast'd  on 
data  largely  from  southern  California  anti  presumed  not  tvpical  of  all  regions  of 

a") 

the  world,  “ 

25.  liattis,  .1.  C.  (1981)  Regional  modification  of  acceleration  attenuation  func¬ 
tions,  Bull.  Seismol.  Sor,  Am.  71:1309-1321. 


Strong  ground  motion  attenuation  in  the  Utah  area  is  not  well  determined. 
Lacking  adequate  data  upon  which  to  base  an  attenuation  function,  we  used  curves 
derived  for  the  California  area.  The  function  parameters  for  the  peak  acceleration, 
velocity,  and  displacement  attenuation  functions  are  given  in  Table  3.  The 


Table  3.  Peak  Ground  Motion  Attenuation  Function  Parameters 
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acceleration  curve  was  derived  bv  Hat t is,  tnd  tin  remaining  two  functions  were 
•;  i  ) 

evaluated  bv  McGuire."  '  "  These  equations  give  peak  ground  motion  value*  for 
rock  uid  stiff  -oil  sites. 

1.3  -I. ih-lic  . il  II. i/. ml  |{e>i il t - 

I  or  .•  jeli  of  the  ’  in  moii,  I  -  .|  •  .  rii ..  . i  jn  Section  peak  ground  tcceleration 
and  lo.itv  eontour  map-  •.  er>  tin  tied  for  I  tab  using  two  different  lifetime 
period-,  10  a  ir  and  Ml  v.  u-,  mil  Mu  il.  nlatiotis  were  earrii  d  out  at  the 
nQ  pert  e nt  . -  onl  idenee  h ■'  .  1  i  >l  non1  *:ei  -. -t !,■  ne 1 ■  ( 1 1  ^  -  0.  1 ),  I  or  the  1 0  - ve ar  lib  f ime, 
tin  ground  motion  levels  correspond  to  the  u.a-vear  return  period  motions:  for  tin 
50-ve.ar  lifetime,  to  4 75 -year  return  period  motions.  Those  contour  maps  arc 
shown  in  Figures  7  through  14. 

Two  facts  concerning  these  maps  should  be  pointed  out.  First,  the  ground 
motion  levels  arc  for  stiff  soil  and  rock  site's;  they  have  not  been  modified  in  any 
wav  for  local  geologic  conditions.  Particularly  in  the  alluvial  basins  of  western 
l  tail,  surface  ground  motions  could  he  expected  to  be  significantly  different  from 
reported  values.  Second,  the  general  contour  trends  are  determined  by  the  areal 
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extent  of  I'aoli  sour ot*  region.  Any  change  in  the  sourer  region  do finitions  would  be 
exported  to  alter  the  ground  motion  predicted  at  any  particular  location.  The 
limitations  of  available  geologic  data  and  the  diffuse  nature  of  seismicitv  in  this 
region  allow  only  gross  outlines  of  source  regions  to  be  made.  This,  in  turn, 
limits  application  of  these  maps  to  anything  more  than  a  regional  overview. 

Variations  in  the  seismic-  hazard  predicted  by  the  two  seismicitv  models  are 
apparent  in  Figures  7  through  14.  This  is  espeeiallv  true  in  western  >  tali  :uul 
along  the  Wasatch  Front  area  and  expected  from  the  different  manner  in  which  each 
model  distributes  the  seismic  activity  in  this  area.  As  would  be  expected  from  the 
seismicitv  models,  the  seismic  hazard  based  on  the  uniform  model  predicts  higher 
levels  of  ground  motion  overall  and  specifically  in  the  Basin  and  Range  Province 
region  of  western  Utah. 

In  addition  to  these  contour  maps,  annual  ground  motion  risk  curves  were 
evaluated  for  four  sites  that  approximately  correspond  to  Cedar  Citv  CHIP'S,  112"  W), 
Wendover  (4F.N,  114'  W),  Salt  Lake  City -Ogden  (41"  X,  112'  W)  and  Manti-I.a  Sal 
National  Forest  (28' N,  110'W),  The  seismic  risk  evaluated  at  these  locations 
show  the  range  of  hazard  in  the  state  uid  demonstrates  the  effects  of  the  model 
variations.  The  annual  risk  curves  for  these  locations  arc  shown  in  Figures  la 
through  22.  The  ground  motion  levels  for  specific  annual  risks  for  each  site  are 
given  in  Table's  4  through  7.  For  purposes  of  comparison,  the  annual  seismic 
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risk  curves  for  Vandenherg  Air  Force  Base,  California  are  shown  in  Figure  22.“ 
Vandenbcrg  Air  Force  Base  is  located  in  a  region  of  California  that  lias  a  moder¬ 
ately  high  si  'smic  hazard  level. 

1.4  ('.oniposit)'  Design  Response  Speelra 

The  spectral  characteristics  of  ground  motion  are  typically  represented  in  the 

form  of  response  spectra.  These  spectra  represent  the  maximum  response  of  a 

simple,  viscous -damped  harmonic  oscillator  over  a  range  of  natural  periods.  For 

specific  levels  of  damping,  methods  hav  e  been  developed  to  estimate  upper  limit 

■>  u 

response  spectra  given  the  predicted  ground  motions  levels  at  some  site.  “  These 


26.  McGuire,  R.  K,  (1U74)  Seismic  Structural  Response  Risk  Analysis  Incorpora¬ 

ting  l’eak  Response  Regressions  on  Earthquake  Magnitude  uid  Distance, 

M.  I.T,  Dept,  of  Civil  Bug.  ,  Research  Report  R7-4-51. 

27.  Battis,  .T.  C.  ( 1 1)7 !))  Seismic  Hazards  Estimation  Study  for  Vandenherg  VFU, 

I  SAF  Geophysics  Laboratory,  AFGL-TR-7!>-0277,  AUA08245H. 

2!i.  Hays,  IV.  W.  ,  Algermisscn,  S.  T. ,  Espinosa,  A.  F.  ,  I’orkins,  I).  M.  ,  and 
Rinehart,  IV.  A.  (1!)75)  Guidelines  for  Developing  Design  Earthquake 
Response  Spectra,  F.  S.  Army  Construction  Eng.  Research  Laboratory, 
Technical  Report  1U-114, 
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Figure  15.  Annual  Seismic  Risk  < 'u e.->  for  1 1  >»•  Cedar  City  Area  l  sinu  l  niform  Seismicity 


Siijure  l<i.  \nnual  Seismic  Hisk  Curves  for  the  Cellar  City  \rea  l  .sin^  Subplate  Margin  Seismicity  Model 


Figure  18.  Annual  Seismic  l{i.->k  Curves  for  the  Uendover  Area  I'sinu  Subplate  Seismicity  Model 


xsm  "Ivdniw 


MB 


J 


xsm  ivhnnv 


DISPLACEMENT  \  \  VELOCITY  \  ACCELERATION 


DISPLACEMENT  \  \  VELOCITY  \  ACCELERATION 


6 

o 


1  ‘  ' _ I _ !■  i  I  i  i  l  i  i _ i  i  i _ ''''i  I _ I _ L 


O 


xsiy  ivnNNv 


o 

o 


o 

o 


u 


T. 

£ 

o 

•H 

E 

V) 

Q) 

CO 


3  U 

■2fJ 

h.  S 


40 


gin  Seismicity  Model 


PEAK  GROUND  MOTION 


ANNUAL  SEISMIC  RISK  CURVES  FOR  VANDENBERG  AFB 


Figure  23.  Annual  Seismic  Ilisk  Curves  for  Pt  Arguello  (Vandenberg  AFB) 
California 


Table  6.  Peak  Ground  Motion  Annual  ltisk  Levels  for  Salt  Lake  City  -  Ogden  Area 


Table  7..  Peak  Ground  Motion  Annual  Hisk  Levels  for  the  Manti-La  Sal  National 
Forest  Area 


are  known  ->s  design  response  spectra.  Newt  lark  et  al  have  developed  one  sot  of 

)  >  j 

commonly  used  amplification  factors  given  in  Table  (>. "  These  values  are  used  to 


Table  fi.  Horizontal  Design  Response  Spectra  Amplification  f  actors 
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modify  the  ground  acceleration  and  displacement  levels  in  order  to  obtain  response 
spectral  levels  at  specific  frequencies.  The  levels  of  critical  damping  correspond 
to  different  foundation  soil  conditions,  with  the  lowest  critical  damping  values  being 
for  hard  rock,  and  the  highest  damping  values  for  soft  soil.  These  spectra  repre¬ 
sent  the  mean  plus  one  standard  deviation  spectral  values  for  any  specified  ground 
motion. 

Composite  response  spectra  have  been  calculated  for  the  100-vear  return 
period  ground  motions  at  each  of  the  four  sites  specified  in  Section  -4.  H.  These 
spectra  are  given  in  figures  24  to  27.  These  curves  correspond  closelv  to  the 
10-vear  lifetime,  !)0  percent  confidence  level  response  spectra.  It  should  be  noted 
that  the  curves  do  not  necessarily  represent  the  predicted  spectra  tor  anv  one 
1  arthquake.  The  peak  ground  motions  at  the  100 -year  return  period  could  be 
generated  bv  different  earthquakes  occurring  at  different  times  and  locations. 

Thus,  it  is  more  accurate  to  view  this  representation  as  estimated  upper  level 
values  for  .uiv  one  event  over  limited  frequency  windows. 


2tl.  New  mark,  N.  M.  ,  Illume,  J.  A.,  and  Kapur,  K.  K.  (1!)7:D  Design  Response 
Spectra  for  Nuclear  Power  Plants,  Am.  Soc.  Civil  frig.  ,  Structural  ring. 
Aleeting,  San  Francisco,  California. 
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Figure  24.  100-Year  Return  Period  Composite  Response  Spectra  for  the 
(  edar  City  Area:  (a)  l  inform  Seismicity  Model;  (b)  Subplate  Margin 
Seismicity  Model 
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FREQUENCY  CHZ) 


Figure  27.  100-Year  Return  Period  Cbmposite  Response  Spertra  for  the 
Manti-T.a  Sal  National  Forest  Area  (both  models  produce  t lie  same  spertra) 


5.  MAXIMUM  ('.RKDIBI.K  GROUND  MOTION 

With  use  of  the  map  of  known  and  suspected  Quaternary  faults  in  1  t.nh  (l  ieure  :t>. 
maximum  credible  ground  motions  were  estimated  for  the  (  tab  region.  As  m.tm 
fault  segments  on  this  map  could  be  incompletely  plotted,  each  fault  segment  was 
assumed  to  be  able  to  support  the  maximum  rredibh  earthquake  for  the  source 
region  in  which  it  is  located.  The  maximum  credible  earthquakes  are  based  on 
the  subplate  margin  seismicity  model  (  Table  2).  Ground  motion  distrinution  is 
based  on  the  00  percent  confidence  level  motions  predicted  hv  the  attenuation  func¬ 
tions  given  in  Table  3.  The  resulting  contour  plots  arc  shown  in  figures  2h  through 
30.  Rased  on  the  available  fault  data,  which  in  m.uiv  areas  is  incomplete,  these 
contour  maps  show  the  largest  ground  motions  expected  al  a  hard  rock  site  in  I  tab 
without  regard  to  the  likelihood  of  occurrence  in  anv  time  period. 


4!! 


Figure  28.  Deterministic  Peak  Acceleration  Contours  for  1  tali  in  Fractions 
of  ';,  (contours  above  0.8  i»  are  not  plotted) 
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(>.  DISCI  SSION  OK  KI.SI  ITS 


The  level  of  uncertainty  in  the  evaluation  of  seismic  hazard  for  any  region  is 
high,  and  particularly  true  for  this  study.  The  source  of  this  uncertainty  includes: 
the  definition  of  the  seismicity  characteristics  throughout  the  region;  the  seismic 
attenuation  properties  of  the  crust;  and  the  relation  between  the  accelerations  gen¬ 
erated  l)v  an  earthquake  and  the  size  of  the  event.  Kite  effect  of  this  uncertainty 
ran  be  seen  by  comparison  of  the  annual  risk  levels  given  in  Tables  4  through  7. 
Variations  greater  than  -100  percent  can  he  found  where  the  only  modification  being 
considered  is  the  distribution  of  seismicity. 

In  the  present  study,  the  definition  of  the  seismic  regionalization  and  tic-  seis¬ 
mic  characteristics  of  each  region,  such  as  maximum  magnitude  and  recurrence 
rates,  are  believed  to  be  the  ma.ior  source  of  error.  As  one  example,  tin-  b-value 
of  the  recurrence  eurvos  is  often  found  to  have  significant  variations  when  derived 
using  distinct  time  windows  of  earthquake  historv  within  a  specified  region,  1  his 
is  a  real  effort  which  derives  from  the  statistical  nonstation aritv  of  the  eai  thou  ike 
proves  s. 

Because  the  length  of  the  historical  record  is  short  when  compared  to  the 
duration  time-scale  of  geologic  processes,  choice  of  the  proper  seismieitv  model 
is  a  difficult  task.  To  some  degree,  the  analysis  done  for  this  report  attempts  to 
handle  those  questions  by  using  two  models  representing  approximations  oi  the 
extremes.  The  subplate  margin  seismicity  model  is  considered  the  better  model  it 
one  projects  historical  seismic  aetivitv  into  the  future.  However,  it  the  historic  nl 
seismicity  does  not  represent  the  long-term  trend,  as  suggested  hv  l’v  ill.  then 
the  uniform  seismicity  could  be  more  accurate. 

l  or  predicting  seismic  hazard  over  the  in  ar-term,  up  o  sever  il  huiidrci 
vears,  the  author  holds  that  the  subplate  margin  model  provides  t  i *>*ti ■  r  picture 
of  future  earthquake  activity.  The  primary  reason  for  this  conclusion  is  tu  u  ns 
the  geologic  processes  occur  over  long  lime  periods,  eh  inge-  ir.  tile  proc  --,  - 
■.'.ill  also  occur  over  relatively  lone  periods. 


7.  CONCl.l  SIONS 

With  the  use  of  two  distinct  seismic  regionalizations,  seismic  h  u-  ird  •  \  1  i- 
tions  were  conducted  for  I  tall,  depicting  the  possible  extremes  in  future  sei-mic 
activity  in  the  1  tali  area.  The  primary  distinction  in  the  models  w  i<  the  tin  nm  ¬ 
ol'  II  as  in  and  llange  seismieitv.  In  one  model,  the  seismic  netivitv  in  tlii-  proven- 


was  uniformly  distributed,  whereas  in  the'  second,  it  was  concentrated  at  the  prov¬ 
ince  boundaries.  On  the  basis  of  these  models,  contour  maps  of  peak  accelera¬ 
tion  and  velocity  were  constructed  for  1'tah  and  annual  risk  curves  evaluated  for 
specific  sites  within  the  state.  In  addition,  deterministic  hazard  estimates  were 
made  based  on  known  and  suspected  Quaternarv  faults  in  rtah. 

It  has  been  concluded  that  for  near-term  hazard  evaluations,  the  subplate  mar¬ 
gin  model  provided  the  most  useful  hazard  projection,  liased  on  tliis  model,  90 

2 

percent  confidence  level  accelerations  in  1  tali  are  less  than  11a  em/sec  for  :uiy 
50-year  period.  The  deterministic  modeling  predicts  that  accelerations  over 
0.  8  g  are  possible  throughout  much  of  western  1‘tah  and,  in  particular,  the  pro¬ 
posed  MX  basing  area.  However,  it  is  likely  that  accelerations  of  this  magnitude 
would  occur  on  the  order  of  less  than  once  in  several  thousand  years  at  any 
specific  site  within  the  state. 
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